Ultracold spin-polarized mixtures of molecules with S'-state atoms: 
Collisional stability and implications for sympathetic cooling 
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The prospects of sympathetic coohng of polar molecules with magnetically co-trapped alkali- 
metal atoms are generally considered poor due to strongly anisotropic atom-molecule interactions 
leading to large spin relaxation rates. Using rigorous quantum scattering calculations based on ab 
initio interaction potentials, we show that inelastic spin relaxation in low-temperature collisions of 
CaH(^E) molecules with Li and Mg atoms occurs at a slow rate despite the strongly anisotropic 
interactions. This unexpected result, which we rationalize using multichannel quantum defect theory, 
opens up the possibility of sympathetic cooling of polar molecules with alkali-metal and alkaline- 
earth atoms in a magnetic trap. 
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The experimental realization of molecular ensembles 
cooled to temperatures below 1 K has opened up a mul- 
titude of fascinating research directions in physics and 
chemistry [ij . The unique properties of ultracold molec- 
ular gases such as their long-range, anisotropic interac- 
tions may be used to implement quantum logic gates 3 , 
create and explore novel phases of quantum matter 3|, 
and study non-equilibrium dynamics and many-body lo- 
calization phenomena Recent experimental and the- 
oretical work has demonstrated the possibility of con- 
trolling chemical reactions in an ultracold gas of KRb 
molecules by applying external electric fields and confin- 
ing the reactants in low dimensions 01, opening up a 
vast new area of research in physical chemistry [ij. 

A variety of experimental techniques has been devel- 
oped for cooling and trapping polar molecules Indi- 
rect cooling techniques such as photo- and magnetoasso- 
ciation [ij assemble ultracold molecules from precursor 
atoms using laser light and oscillating magnetic fields. 
While capable of producing high phase-space density 
gases of ground-state polar molecules, the indirect tech- 
niques are currently limited in scope to E-state alkali- 
metal dimers Alternatively, molecules can be cooled 
directly via thermal contact with cryogenic He buffer gas 
(buffer- gas cooling P, 0|) or molecular beam decelera- 
tion The direct cooling techniques are ideally suited 
for the production of diverse classes of polar molecules 
required for specific applications in quantum simulation 
[H and precision measurements However, these 

techniques tend to produce cold (T > 10 mK) rather than 
ultracold molecules, necessitating an additional stage of 
cooling to reach ultralow temperatures. 

Sympathetic cooling is arguably the most straightfor- 
ward way of cooling molecules below 10 mK by bring- 
ing them in thermal contact with a reservoir of ultracold 
atoms such as the alkali- metal atoms ll-l3- As this 
technique does not rely on specific details of molecular 
structure, it can potentially lead to the production of 



a wide array of ultracold molecules. Molecular cooling 
and trapping experiments typically employ static elec- 
tromagnetic traps, which confine molecules in their low- 
field-seeking Zeeman states Inelastic collisions with 
buffer-gas atoms lead to spin relaxation and trap loss 
[l^ , and a large ratio of elastic to inelastic collision rates 
(7 > 100) is required for efficient cooling [l|. Previous 
theoretical work has shown that spin relaxation in colli- 
sions of alkali- metal atoms with OH [ll|, NH [l3|, and 
NH3 [l3| is extremely fast (7 ^ 1). Based on these re- 
sults, the alkali-metal atoms were judged unsuitable as 
collision partners for sympathetic cooling of molecules, 
and attention turned to the alkaline-earth [l^ , and rare- 
gas [3| atoms, which lack magnetic moments and are 
thus challenging to accumulate in a magnetic trap. Spin- 
polarized nitrogen atoms have favorable collisional prop- 
erties with NH molecules [l^l but their use in sympa- 
thetic cooling experiments is limited by a lack of efficient 
detection techniques. 

In this Letter, we explore the possibility of sympathetic 
cooling of diatomic molecules with one unpaired electron 
by collisions with spin-polarized alkali-metal and alkali- 
earth atoms in a magnetic trap. Using converged quan- 
tum scattering calculations based on ab initio interac- 
tion potentials, we demonstrate that spin relaxation in 
collisions of CaH(^I]) molecules with Li and Mg atoms 
occurs slowly despite the fact that Li-CaH and Mg-CaH 
interactions are extremely strong and anisotropic. To ex- 
plain this surprising result, we map out the dependence 
of the inelastic cross section on the potential anisotropy 
using multichannel quantum defect theory. Our analysis 
shows that spin relaxation in low-temperature collisions 
of molecules can be slow even in systems with strongly 
anisotropic interactions such as Li-CaH, suggesting that 
sympathetic cooling of ^E molecules with alkali- metal 
atoms in a magnetic trap is likely to be successful. 

We begin by outlining our theoretical approach. The 
Hamiltonian of the collision complex formed by a ^E 
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molecule A and an S-state atom B is (fi, = 1) [isl. \li 
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where n is the reduced mass, R is the atom-molecule sep- 
aration, r is the internuclear distance, 9 is the angle be- 
tween the Jacobi vectors R and r, I is the orbital angular 
momentum for the collision, Vgi; (i?, r, 9) is the potential 
energy surface (PES) for the atom-molecule interaction, 
and S = Sa+Sb is the total spin. The Hamiltonian of the 
molecule is given by Ha = BeN'^+jN-SA+2fioBSAz, 
where N is the rotational angular momentum of the 
molecule, B is the magnetic field strength, Saz is the pro- 
jection of molecular spin Sa on the field axis, Be is the 
rotational constant, 7 is the spin-rotation constant (20| . 
and ^0 is the Bohr magneton. The atomic Hamiltonian 
is given by Hb = 2fiQBSBz and the term proportional to 
R~^ represents the magnetic dipole interaction. Both of 
these terms are absent for the Mg atom bearing no mag- 
netic moment. In this work, we are interested in collisions 
of CaH molecules with Li atoms initially in their fully 
spin-polarized Zeeman states Mg^ = Ms^ = 1/2, where 
Msj^ and Mg^ are the projections of Sa and on the 
magnetic field axis. We can therefore neglect the weak 
magnetic dipole and spin-rotation couplings between ba- 
sis states of different S and include only the 5 = 1 
(triplet) Li-CaH PES in scattering calculations. 

The interaction PESs for Li-CaH and Mg-CaH were 
evaluated ab initio jlH using a highly correlated 
open-shell coupled cluster CCSD(T) method and large 
quadruple-^ correlation-consistent basis sets augmented 
by atomic-centered diffuse and i?-centered bond func- 
tions. The CaH fragment was kept frozen at equilibrium 
(re = 3.803ao). To study the dependence of the Li-CaH 
interaction energy on the CaH stretching coordinate r, 
we calculated the PES over the range of r = 3.0— 5.7ao. 
We found that the interaction potential becomes slightly 
deeper with increasing r, but retains its repulsive wall. 
This demonstrates that the chemical reaction Li + CaH 
— >■ LiH + Ca, which is exothermic by ~0.6 eV, is forbid- 
den when both the reactants are spin-polarized. We note 
that the reaction has been observed to occur rapidly at 
cryogenic temperatures with unpolarized reactants [22 1. 

A contour plot of our ab initio PES for Li-CaH is shown 
in Fig. [1}l. The global minimum is located at i? = 5.6 
ao, 9 = 35° and is 7063 cm~^ deep. Two saddle points 
occur in the linear Li-CaH (Li-HCa) configurations at 
R = 10.7ao (7.0ao) with the corresponding well depths 
of —158.5 and —4384 cm~^, respectively. This demon- 
strates that the Li-CaH interaction is not only strong, 
but also strikingly anisotropic. For Mg-CaH, in addition 
to the global minimum of 2110 cm~^ at i? = 6.08ao, 
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Fig. 1. Contour plots of the ab initio PESs for Li-CaH {S — 1) 
(a) and Mg-CaH (b). Energies are in units of cm~^. Collinear 
arrangements of atoms for 6 = Q° and 180° are indicated 
between panels, where Me = Li, Mg. 



9 = 42°, we find a local minimum at the linear Mg-CaH 
arrangement at i? = 7. Gap with the well depth of 1855 
cm~^. As shown in Fig. [1] both the Li-CaH and Mg-CaH 
PESs have similar topology at long range, but differ sig- 
nificantly at small R, with the Li-CaH interaction being 
significantly more anisotropic. 

In order to assess the prospects for sympathetic cooling 
of CaH molecules with Li and Mg atoms in a magnetic 
trap, we solve the quantum collision problem specified by 
the Hamiltonian (jlj numerically using a close-coupling 
(CC) approach in the body- fixed coordinate frame [l^. 
Due to the large anisotropy of Li-CaH and Mg-CaH in- 
teractions, a large number of rotational channels must be 
included in the basis set to obtain converged results. As 
shown in the inset of Fig. [2^, the cross sections oscillate 
dramatically with increasing the size of the rotational 
basis set until convergence is reached at iVmax — 55. 
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Fig. 2. (a) Cross sections for spin relaxation in Li-CaH colli- 
sions calculated as functions of collision energy for B — 0.1 T. 
Full lines and symbols - calculations without the magnetic 
dipole interaction, dashed line - calculations including the 
magnetic dipole interaction, (b) Thermally averaged ratio of 
elastic to inelastic collision rates for Li-CaH as a function of 
temperature. The inset shows the convergence of the inelastic 
cross section at e = 0.05 cm~^ with the maximum number of 
rotational channels included in the basis set (A^max). 

Four total angular momentum states (J ^ 0.5 — 4.5) 
were included in the basis to produce the cross sections 
converged to <10 %. The maximum number of chan- 
nels was 3250 for the total angular momentum projection 
M = 1/2. We note that the same problem reformulated 
in the space-fixed representation |18| with Nj^ax = 55, 
^max = 55 would amount to 234136 channels, making it 
computationally intractable. 

Figure[2Ua) shows the cross sections for spin relaxation 
and elastic scattering in Li-CaH collisions as functions of 
collision energy calculated for a fixed value of the mag- 
netic field {B = 0.1 T). The inelastic cross sections re- 
main small (below 0.5 A) over the entire range of colli- 
sion energies considered except in the vicinity of a shape 
resonance at e ~ 0.05 cm~^. The magnetic dipole in- 
teraction has a dramatic effect on inelastic cross sections 
at low collision energies; the cross sections assume their 
characteristic threshold behavior at a higher e than they 
would in the absence of the interaction. In this regime, 
our calculated spin relaxation rates for Li-CaH are sim- 
ilar in magnitude to those observed for spin-polarized 
alkali-metal atoms At e > 0.001 cm~^, the magnetic 
dipole interaction has little effect on the inelastic cross 
sections, which are dominated instead by the intramolec- 



Fig. 3. (a) Cross sections for spin relaxation in Mg-CaH colli- 
sions calculated as functions of collision energy for B = 0.1 T 
(diamonds) and 1 T (triangles); (b) Thermally averaged ratio 
of elastic to inelastic collision rates for Mg-CaH vs tempera- 
ture for B = 0.1 T (full line) and 1 T (dashed line). 



ular spin-rotation interaction. The Mg-CaH interaction 
is more weakly anisotropic and the cross sections for spin 
relaxation in Mg-CaH collisions shown in Fig. [3^ are 
smaller than for Li-CaH. The ratios of elastic-to-inelastic 
collision rates for Mg-CaH displayed in Fig. [HId are above 
lO'^ at all temperatures, suggesting excellent prospects for 
sympathetic cooling of CaH molecules with Mg atoms. 

Given the extremely large anisotropy of the Li-CaH 
and Mg-CaH interactions, it is remarkable that the in- 
elastic cross sections shown in Figs. [5] and [3] are so small 
in absolute magnitude. In order to gain insight into the 
mechanism of spin relaxation in strongly anisotropic col- 
lision systems, we employ the powerful formalism of mul- 
tichannel quantum defect theory (MQDT) We begin 
by partitioning the full _ftr-matrix into closed and open- 
channel blocks [25| using a minimal basis set consisting 
of two open {\N = 0,Msj^ = ±1/2)) and two closed 
(|iV = l,Msj^ = ±1/2)) channels. In the absence of the 
magnetic dipole interaction, there is no direct coupling 
between the open channels [l3|, and we find K°° = 
and K°'^ = K'^° = dmg{Ki,Ki)^ where Ki quantifies the 
strength of the anisotropic coupling between the TV = 
and = 1 rotational states. The closed-closed block 

of the A'-matrix is given by K'^'^ — ( l^'^ ] , where 
^ \ Kbb J' 

Kbb is the closed-channel level shift [25| and is the 

spin-rotation coupling between the closed channels. We 
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note that is small, and can thus be approximated as 
= —^dv},/de, where V}, is the quantum defect [25| . 
Having defined our model 4x4 if-matrix, we obtain 
the reduced _^-matrix using if'"'^ = K°° + K°''[tanub + 
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where D = (tan Vb + Kbb) 
obtain for the spin relaxation probability 



(2) 

K'l. Using this result, we 
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This expression immediately yields the scaling law 
(Ti/2->-i/2 ~ 7^ [3 3.nd illustrates several important fea- 
tures of the mechanism of collisional spin relaxation in 
molecules. First, the inelastic probability passes through 
a maximum at Ki = [(tanvb + KbbY ~ K'i^l'^ and tends 
to zero as when — > oo. The maximum possi- 

ble value of |S'i/2^_i/2p is given by K'^l\\.ax\Vb + Kb\^ . 
This implies that spin relaxation can be slow even when 
anisotropic interactions are extremely strong {K\ — > oo). 
The physical meaning of this result is as follows. In a 
strongly anisotropic collision system, the interaction po- 
tential mixes the incident iV = collision channel with 
many closed A*" > channels. However, spin relaxation 
cannot occur without the spin-rotation interaction flip- 
ping the electron spin within each iV > manifold. Be- 
cause the spin-rotation interaction is weak, this process 
is inefficient, acting as a "dynamical bottleneck" which is 
qualitatively similar to the role of the transition state in 
abstraction chemical reactions (2^ . 

In the opposite limit of weakly anisotropic interaction 



(A'l ^ 0), Eq. © yields ISi/s-,-! 



1/2 I 



This re- 



sult can also be obtained using third-order distorted- wave 
approximation [2^. Away from Feshbach resonances, the 
lowest non-vanishing contribution to the S'-matrix may 
be written as S'i/2_j._i/2 ^ ^iill^l.-, where Vbi is the 
matrix element of the interaction potential between the 
incident (iV = 0) and closed-channel (iV = 1) distorted 
waves. The cross sections for spin relaxation thus scale 
as V^'\7^/Sg, the result obtained numerically in 

In summary, we have presented a theoretical analysis 
of low-temperature collisions of CaH(^S) molecules with 
Li and Mg atoms based on accurate ah initio interaction 
potentials and a rigorous quantum mechanical approach 
[19| . We have found that the interactions between Li and 
CaH are extremely strong and anisotropic, but the chem- 
ical exchange Li + CaH — > LiH + Ca is not allowed when 
both collision partners are fully spin-polarized. Our cal- 
culations show that inelastic spin relaxation in collisions 
of fully spin-polarized CaH molecules with Li and Mg 
atoms occurs slowly, with elastic-to-inelastic ratios in ex- 
cess of 10^ over a wide range of temperatures from 10 



/iK to 10 mK. We developed an analytic model, which 
predicts that spin relaxation cross sections scale as K\ 
with the interaction anisotropy for small K\ , but decrease 
to zero as K\ — >■ oo. This non- linear scaling should be 
applicable to all collision-induced relaxation processes in- 
volving molecules at temperatures below the first ro- 
tational excitation threshold [kT < 2Be). 

Our results may open up novel avenues of research with 
cold and ultracold molecules. First, they demonstrate 
that ultracold spin-polarized mixtures of molecules 
with alkali-metal atoms are stable against collisional spin 
relaxation, strongly suggesting that sympathetic cooling 
of polar CaH(^E) molecules with Li and Mg atoms will 
be successful. It would be interesting to explore the 
prospects of sympathetic cooling of laser-cooled SrF(^S) 
and electric-field-decelerated YbF(2E) Second, 
our analysis indicates that collisions of molecules with 
each other are likely to be predominantly elastic, and 
hence evaporative cooling of molecules in a magnetic 
trap may be feasible. Finally, our work opens up the 
possibility of sympathetic cooling of polyatomic spin- 1/2 
radicals, which can be magnetically trapped at milli- 
Kelvin temperatures via buffer-gas cooling (ist . 
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